A dynamical law for slow crack growth in polycarbonate films 
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We study experimentally the slow growth of a single crack in polycarbonate films submitted to 
uniaxial and constant imposed stress. For this visco-plastic material, we uncover a dynamical law 
that describes the dependence of the instantaneous crack velocity with experimental parameters. 
The law involves a Dugdale-Barenblatt static description of crack tip plastic zones associated to an 
Eyring's law and an empirical dependence with the crack length that may come from a residual 
elastic field. 

PACS numbers: 62.20.Mk, 62.20.Fe, 46.35.+Z 



Stressed solids commonly break apart once a critical 
stress threshold is reached. However, many experiments 
0) S S BL HI show that a given solid submitted to a sub- 
critical stress breaks after a certain amount of time. The- 
refore, understanding the mechanisms of subcritical rup- 
ture of solids has become an important goal of fracture 
physics in order to improve the resistance of structures to 
delayed failure that may have catastrophic consequences. 
According to reported experimental works [H, [3] , the de- 
pendence of the rupture time with applied stress a can 
be described in many kinds of materials (polymers, me- 
tal alloys, semi-conductors, rocks...) by an Arrhenius law 
with an energy barrier decreasing linearly with a. This 
proposed universality is disturbing since these materials 
have micro-structures and rheological properties very dif- 
ferent from one another, and the rupture dynamics is 
certainly expected to be dependent on those properties. 
To lift this paradox, one must go beyond characteriza- 
tion of global properties such as rupture time and ins- 
tead study experimentally the full time-resolved rupture 
dynamics, from the stress application to the final break- 
down of the sample. A convenient system to start with is 
a two-dimensional solid with a single macroscopic initial 
crack submitted to a uniaxial constant load. 

In this context, recent experimental studies H have 
shown that subcritical crack growth in paper sheets 
can be successfully described by a thermally activated 
mechanism inspired from previous theoretical works in 
elastic brittle media @, 0] ■ Experimental study of slow 
crack growth in a visco-plastic material under stress is 
a very active topic @|. General theoretical frameworks 
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ll| have been proposed to predict the dependence 



of the crack growth velocity with experimental parame- 
ters using characteristic material time-response functions 
such as its compliance. However, these models involve 
complex integro-differential equations which are hardly 
tractable in practical situations where visco-plastic ef- 
fects are strong. Consequently, the experimental time 
evolution of the instantaneous crack growth dynamics 
can not be captured easily by current models. 

In order to provide more experimental insight in our 



Fig. 1: Image of a crack in a polycarbonate film with its 
macroscopic plastic zone at each tip. 



understanding of visco-plastic effects during slow crack 
growth, we have performed an experimental study of 
the slow growth of a single crack in amorphous poly- 
mer films made of polycarbonate which is a highly non- 
brittle visco-plastic material. The experiments consist in 
the growth of a single linear crack in a polycarbonate film 
submitted to uniaxial and constant imposed force. The 
polycarbonate films used are Bayer Makrofol® DE and 
have the properties of bulk material. Before each expe- 
riment, a crack of length £i (from 0.5 to 3cm) is initiated 
at the center of the polycarbonate sample (height 21cm, 
length 24cm). Then, a constant force F is applied to the 
film perpendicularly to the crack direction, so that we get 
a mode 1 crack opening configuration. Using a camera, 
we follow the growth of the crack length I under constant 
applied stress a — F/eH (e is the film thickness and H 
the sample height) until the total rupture of the sample. 
The applied stress a is chosen such that crack growth 
is slow, i.e. smaller than a critical one o~ c , above which 
crack propagation occurs in a few seconds. More details 
about the experimental set-up can be found in [l2| ■ 

In each crack growth experiment, during the loading 
phase of the film, a macroscopic flame-shaped plastic 
zone appears at each tip of the crack [lj| and grows with 
the applied stress (cf. Fig. Q] where is defined the plastic 
zone length from tip to tip £ pz ). In the late loading stage, 
the crack may also start to grow at a time that appears 
to be statistical. It is probably a consequence of the dis- 
persion in the local toughness of the material or in the 
initial crack tip shape. Consequently, the real experimen- 
tal initial condition, obtained when the constant stress a 
is reached, is not exactly £ = ij,. Depending on the mo- 
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ment when the crack starts to grow during the loading 
phase, the true initial condition of the creep experiment 
will be a couple of value for the crack and plastic zone 
length (Fig. [1]) : (£*, £* z ). Finally, during the imposed 
stress stage, the plastic zones and the crack are both 
growing until the final breakdown of the sample in a way 
that the crack never catches up the plastic zone tip. In- 
side the plastic zone, the film is subjected to a thinning 
which brings its thickness from 125/im to about 75±5/im. 
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Fig. 2: Time as a function of both the crack and process 
zone lengths for an imposed stress experiment {ti = 1.5cm, 
F — 900N). We indicate the position of the inflexion point t x , 
l x of the crack growth curve. 



Typical growth curves of the fracture and plastic zone 
are shown in Fig. [2] Both curves show a quite similar 
smooth shape. This regular shape lets us think that the 
crack growth in polycarbonate films is a deterministic 
phenomenon. However, for identical experimental condi- 
tions, we notice a large dispersion of the rupture times 
and more generally of the crack growth dynamics. There 
is actually up to a factor five between the rupture time 
of the fastest and slowest experiments. We suggest that 
the explanation for this statistics in the crack growth dy- 
namics does not come from the growth mechanism itself, 
but is a consequence of the dispersion in the effective 
initial conditions at the beginning of the constant stress 
phase of the experiment (£*, ^* z ). These initial condi- 
tions are clearly statistical and hardly controllable in our 
experiment. They are dependent on the moment when 
the crack starts growing during the loading stage of the 
sample and they determine all the rest of the experiment. 

In Fig. [3j we show the evolution of the average rupture 
time (T r ) (averaged over at least ten experiments) as a 
function of the applied stress for a series of experiments 
performed at li — 1.5cm. We see a linear dependence of 
log(T r ) with the applied stress that corresponds well to 
an exponential description of the rupture time as pro- 
posed by Zhurkov [![. The linear fit of the data is of 
quite good quality and suggests that (T r ) — T a e~ arT . In 
Zhurkov's approach, the stress dependence of (T r ) is in- 
terpreted as an Eyring's law [3| with a = V/ksT where 
V is assumed to be a characteristic volume of the mate- 




FlG. 3: Natural logarithm of the average rupture time as a 
function of the applied stress for a series of experiments per- 
formed for ti = 1.5cm. 



rial. However, in our experiments, the parameter V can 
not be a constant since a different initial crack length £i 
gives a completely different rupture time for the same ap- 
plied stress. Thus, the external applied stress a can not 
be the single control parameter of the rupture dynamics. 
Then, it is clear that Zhurkov's description needs to be 
improved to take into account the specific geometry of 
the problem. In particular, the stress a y holding in the 
plastic zone close to the crack tips most probably parti- 
cipates in the dynamical processes leading to the crack 
growth. 
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The Dugdale-Barenblatt cohesive zone model 
is a good and simple mean to estimate the stress a y . This 
quantity appears intuitively as a possible control parame- 
ter for the crack dynamics just like the stress intensity 
factor is for brittle materials. Dugdale-Barenblatt model 
predicts : 
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This plastic stress a y can be computed at each moment 
using Eq. JT]) with the instantaneous values of a, I and 
£ pz . To account for the global dynamics during an ex- 
periment, we compute the time- averaged growth velo- 
city on the whole experiment v and compare it to the 
time-averaged plastic stress ~dy (see Fig. [4]). Each point 
of this Figure represents the mean behaviour over an ex- 
periment. The data are compatible with a linear law that 
predicts an exponential dependence of the average growth 
velocity with the mean stress in the plastic zone : 



v e" 



and v = 7.810- 21 m.s" 



(2) 



It is striking that the prefactor of the stress in the 
exponential curve (cf. Eq. ((2])) is close quantitatively to 
the one obtained in the Eyring's law for the polycarbo- 
nate creep [l2| • Both prefactors probably correspond to a 
unique material constant V/ksT. Thus, we can conclude 
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Fig. 4: Natural logarithm of the average crack growth velo- 
city v as a function of the average plastic stress during the 
growth. Each point represents the average dynamical beha- 
vior during an experiment. Experimental conditions are va- 
rious (£i = 1.5,2,3cm and 2.9 < a < 3.8 10 7 N.m~ 2 ). Each 
experimental condition corresponds to different symbols. 
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Fig. 5: Natural logarithm of the instantaneous crack growth 
velocity as a function of (a) the Dugdale-Barenblatt stress, (b) 
the corrected Dugdale-Barenblatt stress cr™" according to Eq. 
@ for height experiments performed with various experimen- 
tal conditions (£, = 1.5, 2, 3cm and 2.9 < a < 3.8 10 7 N.m" 2 ). 
In Fig. (b), the black line is the result of a linear data fit. 



that the Eyring's law plays a central role in the mecha- 
nisms of crack growth in polycarbonate films. 

We now go beyond a simple analysis of the average 
growth dynamics by looking at the dependence of the 
crack velocity with the stress in the plastic zone at each 
time during the crack growth. We plot in Fig. OUa) the 
instantaneous crack velocity v = d£/dt as a function of 
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Fig. 6: (a) k£ x for various experimental conditions (£i = 
1.5, 2, 3cm and 2.9 < a < 3.8 10 7 N.m~ 2 ) as a function of the 
applied stress a. (b) Critical rupture stress a c as a function 
of the initial crack length li. The dotted lines are a linear fit 
of the data. 



the instantaneous value of the Dugdale-Barenblatt stress 
a y for height experiments performed with various expe- 
rimental conditions. Here, the description of the instan- 
taneous velocity by an exponential law fails, especially 
when the crack length becomes larger than £ x at which 
the minimum crack velocity is reached (£ x is also the in- 
flexion point of the growth curve in Fig. [21). In fact, the 
Eyring's law given by Eq. @ describes well the beha- 
viour only when £ ~ l x . We discovered that introducing 
a correction to a y linear with the crack length £ allows us 
to collapse the experimental data on a straight line (cf. 
Fig-IHb))- This correction can be written as : 
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+ n(£-£ x ). 



(3) 



For each experiment, we determine the value n = (3.4 ± 
0.6) 10 8 N.m~ 3 . The dispersion of k values seems to be 
statistical as no systematic dependence with a or £i could 
be found. This rescaling means that the crack growth 
velocity seems to follow : 



d£ 

d-t =Voe ' 



(4) 



The collapse of the data for various experimental condi- 
tions means that Vq can be considered as a constant. 

In Eq. ([3]), the crack length at the inflexion point in 
the growth curve plays a particular role. It turns out that 
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its value depends on the experimental conditions. This 
can be seen in Fig. EJ^a) where the product k£ x , with 
K = 3.4 10 8 N.m -3 , is plotted as a function of the applied 
stress a. Remarkably, the dependence of k£ x with a is 
well approximated by a linear relation : n£ x = a x — er, 
where a x = 4.2 10 7 N.m~ 2 . 

A way to clarify the meaning of this relation is to look 
at the dependence of the critical stress a c needed to break 
instantaneously a sample with a crack of initial length £{. 
In brittle materials, we would expect this critical stress 
to decrease in 1 /y£~i since the rupture criterion is reached 
when the initial stress intensity factor equals the tough- 
ness of the material K c [|| : a c ^ir£i/2 — K c . For an 
amorphous visco-plastic material such as polycarbonate, 
we do not get the same functional dependence. Indeed, 
as we can see in Fig. &b) , the relation between a c and £i 
can be approximated by a linear equation : (3£i = (t s — <j Ci 
where a s = 4.0710 7 N. m - 2 and j3 = 3.57 10 8 N. m - 3 . We 
note that k ~ /3 and a x ~ a s and will consider these 
quantities to be the same material constants. So, we find 
that the quantity S(cr c , £{) = a c + n£i may play a similar 
role than the initial stress intensity factor in brittle ma- 
terials. Furthermore, it allows us to interpret the value 
of the crack length at the inflexion point as defined by a 
characteristic value of the quantity S(cr, £ x ) — a x ~ o~ s 
that corresponds to an intrinsic property of polycarbo- 
nate. Indeed, a s corresponds to the rupture threshold a c 
in the limit when there is no initial crack. 

According to the previous analysis of the instantaneous 
crack velocity, crack growth in polycarbonate films ap- 
pears to be ruled, during an experiment, by an Eyring's 
law (cf. Eq (J)) with : 
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This effective stress cr™ rr 



is composed of the Dugdale- 
Barenblatt estimation of the crack tip plastic zone stress 
Oy , a linear dependence with the crack length k £ and the 
applied stress at the borders of the sample a. Note that 
in Eq. ^ appears naturally a volume V ~ 2.8 10~ 27 m 3 
close to the one used to describe the simple creep 
flow of polycarbonate (3.1 10 _27 N.m~ 2 ) as well as the 
growth of a necking instability in polycarbonate films 
(3.0 10~ 27 m 3 ). This observation reinforces the idea that 
the Eyring's law for crack growth is truly a consequence 
of the creep behavior of polycarbonate. 

In Eq. ([5]), the viscous relaxation is taken into ac- 
count by the experimentally measured evolution of the 
ratio £/£ pz as the crack grows. Indeed, if this ratio was 
constant, the stress in the plastic zone would also be 
constant and the velocity would increase monotonously 
due to the linear term in crack length. In that case, the 
behavior would actually be qualitatively the same as the 



one for crack growth in brittle facture @. To predict fully 
the viscous dynamics of the crack, we need a second equa- 
tion that will prescribe £ pz : 



d£» 



dt 



f(£ pz ,£,£,a,...) 



(6) 



An original theoretical approach recently developed 
by Bouchbinder [13] in extension to the Shear- 
Transformation-Zone Theory proposed by Falk and Lan- 
ger [l8| is certainly useful for deriving an equation of the 
plastic zone velocity (cf. Eq. ([6])). Additionally, numerical 
simulations that can reproduce the complex visco-plastic 
behavior of polycarbonate may help in goin g fu rther in 
the interpretation of our experimental results lil 20L 21 1 . 



[9 
[10 

[11 

[12. 

[is; 

[14 

[15. 
[16 
[17 

[is; 

[19 

[20 
[21 



784 



3565 



* Electronic address: Loic.Vanel@ens-lyon.fr 
S. N. Zhurkov, Int. J. Frac. Mech. 1, 311 (1965). 

F. Bueche, J. Appl. Phys. 26, 1133 (1955); 28, 
(1957). 

S. S. Brenner, J. Appl. Phys. 33, 33 (1962). 
L. Pauchard and J. Meunier, Phys. Rev. Lett. 70, 
(1993). 

S. Santucci, L. Vanel and S. Ciliberto, Phys. Rev. Lett. 
93, 095505 (2004) ; S. Santucci et al., Europhys. Lett. 
74(4), 595 (2006). 

C. Hsieh and R. Thomson, J. Appl. Phys. 44, 2051 
(1973) ; M. Marder, Phys. Rev. E 54, 3442 (1996). 

S. Santucci et al., Europhys. Lett. 62, 320 (2003). 
N. Haddaoui, A. Chudnovsky and A. Moet, Polymer 27, 
1377 (1986) ; M. T. Takemory, and R. P. Kambour, J. Ma- 
ter. Sci. 16, 1108 (1981) ; S. Hashemi and J. G. Williams, 
Plast. Rubber Compos. 29, 294 (2000). 
R. A. Schapery, Int. J. Fracture 11, 141 (1975) ; 11, 369 
(1975) ; 11, 549 (1975). 

A. A. Kaminskii, Soviet Appl. Mech. 15, 1078 (1979) ; 
Int. Appl. Mech. 40, 829 (2004). 

A. Chudnovsky and Y. Shulkin, Int. J. Fracture 97, 83 
(1999). 

P.-P. Cortet, L. Vanel and S. Ciliberto, Europhys. Lett. 
71(2), 242 (2005) ; in preparation for Eur. Phys. J. B. 

A. M. Donald and E. J. Kramer, J. Mater. Sci. 16, 2967 
(1981). 

H. Eyring, J. Chem. Phys. 4, 283 (1936). 

D. S. Dugdale, J. Mech. Phys. Solids 8, 100 (1960). 

G. I. Barenblatt, Adv. Appl. Mech. 7, 55 (1962). 

E. Bouchbinder, J. S. Langer and I. Procaccia, Phys. Rev. 
E 75, 036107 (2007) ; 036108 (2007). 
M. L. Falk, J. S. Langer and L. Pechenik, Phys. Rev. E 
70, 011507 (2004). 

M. C. Boyce, E. L. Montagut and A. S. Argon, Polym. 
Eng. Sci. 32, 1073 (1992). 

B. P. Gearing and L. Anand, Int. J. Solids Struct. 41, 
827 (2004). 

N. Saad-Gouider, R. Estevez, C. Olagnon and R. Seguela, 
Eng. Fract. Mech. 73, 2503 (2006). 



